Attosecond pulse trains as multi-color coherent control 
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We present a general description of the interaction between multi-color laser pulses and atoms and 
molecules, focusing on the experimentally relevant example of infrared (IR) pulses overlapped with 
attosecond pulse trains (APTs). This formulation reveals explicitly and analytically the role of the 
delay between the IR pulse and APT as a coherent control parameter. Our formulation also shows 
the nearly equivalent roles of the delay and the carrier-envelope phase in controlling the interference 
between different multiphoton pathways. We illustrate these points by investigating the single 
ionization of He and introduce dressed adiabatic hyperspherical potentials to aid the discussion. 
We confirm the predictions with a full-dimensional, two-electron solution of the time-dependent 
Schrodinger equation. 

PACS numbers: 32.80.Rm,42.50.Hz 



At the heart of atomic and molecular physics is the 
motion of bound electrons. Since the relevant time 
scale for their motion is on the order of atomic units 
(1 a.u.«24 as), detailed studies of electron dynamics re- 
quire forces that act on these time scales. With such 
forces, an avenue for control over the electron dynam- 
ics is also created. These abilities to study and control 
are the central goals of the nascent field of laser-based 
attosecond physics [3, Q . 

Attosecond pulses are generated experimentally via 
high-harmonic generation, either as a single isolated 
pulse Q or as an attosecond pulse train (APT) 
The latter, which we will focus on in this Letter, con- 
tains a large range of harmonics of the fundamental in- 
frared (IR) laser (usually a ThSapphire at approximately 
800 nm) . With selective filtering, a specific range of har- 
monics — typically in the extreme ultraviolet range — 
can be isolated and used to illuminate the target. Tem- 
porally, these harmonics combine to give a periodic se- 
ries of pulses, each of which is an attosecond pulse, thus 
the label APT. When overlapped with an IR pulse, an 
APT permits control through the ability to vary the de- 
lay between the APT and IR field which can be done 
on a tens-of-attoseconds scale. Experimental control of 
the ionization probability using such a method has been 
realized E|, Q , as has control over the angular distribu- 
tion @, sT 

Electron dynamics in the combined IR+APT field can, 
of course, be described by solving the time-dependent 
Schrodinger equation [9(. To try to gain some physical 
insight, though, both a time-dependent, electric field pic- 
ture 5|. llOfl and a time-independent, photon based pic- 
ture [111 . Il2j | have been applied. An approach akin to 
the "simple man's approximation" has also been invoked 
to describe the dependence on delay using wave packet 
replicas 

While these attempts to understand the dynamics in 
IR+APT laser fields have some interpretive power, they 
are not general, exact, nor do they provide a simple 



route to predicting the outcome of a particular exper- 
iment. The numerical solution of the time-dependent 
Schrodinger equation is, in principle, general and exact, 
but is limited in practice to simple systems and cannot 
be said to provide a simple guide to predictions. Thus, 
a rigorous description of the problem is desired that can 
also provide a conceptual framework within which pre- 
dictions can be made simply. 

In this Letter, we present just such a theoretical frame- 
work. This formulation is closely related to one we have 
developed [3 Ell to understand the effects of the carrier- 
envelope phase (CEP) for few-cycle pulses and, in fact, 
demonstrates that both the CEP and the IR+APT delay 
provide essentially identical control over a system. This 
conclusion follows because both of these parameters can 
be thought of as controlling the relative phase between 
different multiphoton processes that interfere to produce 
a given physical observable. In fact, attosecond pulses 
can be characterized by taking advantage of such inter- 
ferences [3, 0, E3] • These systems are realizations of the 
canonical coherent control scenario proposed by Brumer 
and Shapiro Moreover, because an APT is a syn- 

thesis of harmonics, our analysis applies equally well to, 
and provides a connection to, the well-studied two-color 
experiments on controlling molecules HI, 3 IS, 21 1 . 



Rather than presenting the most general formulation 
of our theory — which, while conceptually simple, is 
notationally cumbersome — we illustrate it by apply- 
ing it to the single ionization of helium by combined 
IR+APT fields 5]. We confirm our analysis with so- 
lutions of the full six-dimensional, two-electron, time- 
dependent Schrodinger equation. In the process, we thus 
account for the modulation in ion yield with delay be- 
tween IR and APT recently observed in Ref. 5] both 
analytically and numerically. 

Our analytic treatment begins with the time- 
dependent Schrodinger equation in the dipole approxi- 
mation using the length gauge (atomic units are used 
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throughout), 



i-Q*(t)= H -£(t)-d 



(1) 



where Hq is the field-free Hamiltonian and d is the dipole 
operator. We take into account the fundamental IR field 
with frequency u> and two harmonics of order n\ and rii: 



£(t) = £\{t) cos (u>t + <pi) 
+ £ ni (t) cos (ni [w{t - 
+ £ n2 (t) cos (n 2 [w(t- 



(2) 



The latter two terms make up the APT in our example. 
More harmonics can be added to shorten the individual 
pulses in the APT without any conceptual difficulty, but 
three colors is the minimum needed to describe the most 
significant effects observed in an IR+APT experiment 
such as Ref. Q and simplifies the discussion. Note that 
for our analysis the harmonics' pulse envelopes £ ni (t), 
which represent the envelope for the whole pulse train, 
must be independent of the IR-APT delay t but are oth- 
erwise arbitrary. Much like the derivation of CEP effects 
in [l5|], we note from Eq. |T]) that 'J? is periodic in r as 
well as in the IR CEP ipi. We can thus expand \P in both 
using a discrete Fourier transform, 



-Hip 



(3) 



where ip m i (t) does not depend on either ror^i. Substi- 
tuting this W into Eq. ([1]) and equating coefficients of the 
linearly independent exponential functions, we obtain 

d 1 - - 

1 



-l-\-n2 ^ ^ 2 4 ) m — n2,l—n2 



(4) 



In these equations, we have made the additional substi- 
tution i() m i(t) — e~ lluJt (f> m i(t) to emphasize the connec- 
tion with the Floquet representation: in the limit that 
the pulse lengths go to infinity, Eq. ([3]) reduces to many- 
mode Floquet theory [22j . We stress that Eq. (g]) is exact 
even for few-cycle pulses. We could further remove the 
harmonic phases (p rii with Fourier transforms, but that 
will not be necessary for the current discussion. 

We have rigorously obtained a representation in which 
IR+APT experiments can be understood. Further, 
Eq. © shows that the CEP and IR-APT delay pro- 
vide similar control over the dynamics of a system. The 
physical interpretation of this control as interference of 
different photon pathways stems from the facts that the 
components (j) m i oscillate with time-dependence e lluJt and 
correlate to the photon components in the Floquet limit. 
Since we consider photons with different frequencies, the 
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FIG. 1: (Color online) Adiabatic hyperspherical potential 
curves for He. Shown are L = 0, 1, 2 and the region of double 
excitation. 



index I represents only the net energy lu) of the photons 
involved, and can be written as l=h + ini n i+in 2 n 2 with 
li the number of photons of each color. 

Specializing the above discussion to helium determines 
Hq and d in Eq. (j4|). We treat the full three-dimensional 
motion of both electrons using the adiabatic hyperspher- 
ical representation, which has proven especially useful in 
the past for understanding the dynamics of He 2^, 2^ ■ 
More specifically, we expand each component <f> m i as 



= Y J G ™lLv{R 1 t)<S> Lv {R-tt) 



(5) 



where R is the hyperradius and fl the five hyperangles. 
The channel functions <&l„ solve the field-free fixed- R 
adiabatic equation H^^^—Uhv^Lv for a given total or- 
bital angular momentum L and thus form a complete set 
at each R indexed by v. We use Delves' hyperangle and 
treat the system in the body frame, expanding the Eulcr- 
angle dependence of <& Lu on Wigner D-functions (see [24[ 
for details). The adiabatic equation thus reduces to cou- 
pled two-dimensional partial differential equations that 
we solve using b-splines 25|. 

Upon substituting <p m i from (|5|) into Eq. ([4]) and pro- 
jecting out the ^Lv we find coupled, time-dependent 
equations for G m iLv(R,i) describing laser-driven motion 
on the adiabatic potentials Ulw Since the experiments 
in [B[ considered only single ionization at energies well 
below doubly-excited resonances, we can to a good ap- 
proximation consider only the lowest singlet adiabatic 
potential v—1 for each L. These potentials support the 
singly-excited states IsnL 1 L, and several examples are 
shown in Fig.[T] With this truncation of the v expansion, 
we find a ground state energy of -2.895 a.u. 

One advantage of using the adiabatic hyperspherical 
representation in the Floquet-like equations ^ is that 
the resulting dressed potentials — which include the ef- 
fects of both Coulomb interactions and the laser field 
— give insight into the electronic dynamics in much the 
same way that Floquet potentials have for the nuclear dy- 
namics of molecules. Applied to HJ, Floquet potentials 
have provided natural explanations of bond-softening, 
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above-threshold dissociation, vibrational trapping, and 
zero-photon dissociation jig], 27 1. Previous generaliza- 



tions of the Floquet picture for Hj have even helped 
uncover unexpected phenomena such as above-threshold 
Coulomb explosion [28[. It would be interesting to ex- 
plore the analogues of these mechanisms for He in the 
present representation, but we will focus here on their 
application to the IR+APT case. 

Figure [2] shows the diabatic dressed potentials assum- 
ing lj corresponds to 800 nm and the harmonics are 
n! = 13 and ri2=15 (the two most intense components of 
the APT in [5j by a factor of at least 5). These curves 
are generated by first noting that the initial state of the 
system, Is 2 1 S e , has (m,Z)=(0,Q). Then, per Eq. g]), 
this state couples to the 1 P° channel with one IR photon 
(0,±1) or with one harmonic photon (±nj,±nj). The 1 P° 
potential in each case is added to the figure shifted by lui. 
Each of these new curves can couple to either 1 S e or 1 D e 
symmetries according to dipole selection rules, shifted in 
energy by an appropriate amount. Repeating this pro- 
cedure, all possible multiphoton pathways combining IR 
and harmonic photons to any particular final state can be 
generated as shown in Fig. [21 The energy in such a figure 
can be thought of as being approximately conserved, so 
that the dynamics reduces to the well-understood prob- 
lem of multichannel dynamics. This picture, for example, 
maps ionization onto predissociation of molecules in the 
Born-Oppcnheimer approximation or, equivalently, onto 
atomic autoionization in the field-free adiabatic hyper- 
spherical representation. More specifically, transitions 
are most likely to occur where potentials cross, and cross- 
ings of curves differing by the fewest photons will domi- 
nate. In fact, diagonalizing the right-hand-side of Eq. (|U) 
at fixed R will produce adiabatic Floquet-like potentials 
in which these crossings become avoided, thus incorpo- 
rating the strength of the field in the potentials them- 
selves. The remaining time dependence in Si modulates 
the strength of the coupling — or the gap at the avoided 
crossings — as a function of time. 

For the parameters above, single ionization of He re- 
quires the absorption of at least 16w worth of energy, 
i.e. 16 IR photons, 1 13th harmonic plus 3 IR pho- 
tons, or 1 15th harmonic plus 1 IR photons. Since 
it requires the fewest photons, we expect the latter 
scenario to be most likely, leaving the system with 
(to, /)= (15,16) and L—0 or 2 based on dipole selection 
rules. All together, we would write these pathways as 
(0,0) 1 5 e -»(15,15) 1 P°^(15,16) 1 S' e or 1 D e . Following the 
dressed potentials in Fig. [2] suggests that only L=2 will 
be populated substantially, however, since it crosses the 
(15,15) 1 P° curve and L—0 does not. The pathway be- 
ginning with the 13th harmonic would mainly proceed 
via (0,0) 1 S ,e ^(13,13) 1 F°^(13,16) 1 D e by similar argu- 
ments. The above two paths are highlighted in green 
in Fig. [2] Another set of likely pathways (light blue in 
Fig. [2D ends on 1 G e : (0, 0) 1 S e -Kl5,15) 1 P°->(15,16) 1 G e 
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FIG. 2: (Color online) Diabatic dressed potentials for He in a 
three color field (u, 13a;, and 15a;). Only the most significant 
curves and crossing are shown, where a circle indicates a single 
large-photon crossing; a square a 1-IR photon crossing; and 
a triangle, a 3-IR photon crossing. The dashed line indicates 
the energy of the system, which is approximately conserved 
in this representation. 



and (0,0) 1 5 e ^(13,13) 1 P°^(13,16) 1 G e . Other outcomes 
are certainly possible and will occur with some proba- 
bility. The goal of this kind of analysis, though, is to 
try to identify the main routes to a given final state. 
Since there exists more than one pathway to this lowest 
above-threshold ionization (ATI) peak, interference will 
occur. Similar considerations show that there are gen- 
erally multiple paths to ionization for each higher ATI 
peak Z=17,18,19,. . . opening up the possibility of control. 

In this manner, our formulation of the problem lets us 
learn a considerable amount about controlling this inter- 
ference without solving the time-dependent problem. To 
see explicitly how these pathways interfere, we can derive 
an expression for any physical observable using Eq. ([3]). 
In particular, for the He IR+APT photoelectron spec- 
trum discussed above, we obtain 



dE 
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pi(jTt— m')o;Tgi(/ — l)<pi 



x {EL\F mlL ){EL\F m , VL )\ (6) 



where \EL) is an energy-normalized continuum state 
with asymptotic kinetic energy E and F m iL=e~ lluJt G m iL- 
All of the t and ip± dependence is displayed analytically 
since the amplitudes (EL\F m ix) are independent of these 
parameters, determining only which pathways will inter- 
fere. Equation ([6]) shows that the delay dependence is 
determined only by the difference in to, essentially the 
harmonic order, between two pathways. For the lowest 
ATI peak discussed above, the final states of two of the 
pathways identified were (13, lG) 1 !? 6 and (15, ld>) 1 D e for 
the 13th and 15th harmonics, respectively. Since these 
have the same final I (net energy absorbed from the field), 
they will contribute at the same final energy and inter- 
fere, resulting in a sinusoidal dependence on cut with pe- 
riod 7r since |to-to'|=2. This 7T periodicity in the ion- 
ization yield is exactly what was observed in both the 
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FIG. 3: (Color online) (a) Calculated photoelectron spectrum 
for an 800 nm IR pulse with h = 5.0 x 10 12 W/cm 2 and 
7i3 = Ji5 = Clii. (b) Ionization probability as a function of 
delay for L—0, 2, and 4. 

experiment and numerical results presented in [3, 0| . In 
general, the periodicity in lot will be given by 2ir /\m— m'\. 
Thus, if even harmonics had also been included, the pe- 
riodicity would have been 2ir since \m-m'\ would equal 
1. To see the 2tt periodicity, however, one would have 
to measure a quantity that - unlike the energy spectrum 
© - takes into accout interference between different final 
Ls, such as the momentum distribution. 

To confirm our analysis, we carried out a full- 
dimensional numerical solution of the time-dependent 
Schrodinger equation (TTJ using the adiabatic hypersphcr- 
ical representation. Figure [3Ja) shows the resulting pho- 
toelectron spectrum, again forparameters chosen to be 
similar to the experiment in [5j. As expected, the cal- 
culations show ATI peaks separated by to. The largest 
peak is the lowest one (indicated by the arrow), which 
was discussed above, and the predicted it period of its 
modulation is clearly seen. A closer look at the 16w peak 
reveals that it is split in energy with the lower-energy 
peak primarily L=2 and the higher-energy peak L=4 [see 
Fig- EKb)] . These details, however, only serve to bolster 
the value of the dressed potentials we introduced in Fig. [3] 
since they allowed us to predict that the 16w peak should 
be dominated by L—2 and 4. The time-dependent cal- 
culations and the dressed potential picture complement 
each other, providing both quantitative and qualitative 
understanding of these processes. 

This discussion hardly exhausts the insight that can 
be gleaned from Eq. Another example is the fact 

that the delay r between IR and APT will not generate 
a left-right asymmetry along the laser polarization since 
all of the Ls corresponding to a given ATI peak are either 
even or odd (with only odd harmonics in the APT). An 
exception to this statement occurs for broad bandwidth 
pulses since the interference between even and odd Ls 
necessary for asymmetry can be produced in the overlap 
of neighboring ATI peaks. Of course, the CEP will also 
contribute to the asymmetry in this case. 

In this Letter, we have presented a general, exact rep- 
resentation for the interaction of multicolor laser pulses 



with atoms and molecules upon which a relatively simple 
interpretive, predictive picture can be built. In the pro- 
cess, we introduced dressed potentials for He — again 
based on a rigorous derivation — that confer many of 
the benefits such potentials have provided for molecules 
on the atomic problem. We have also shown that con- 
trol via the CEP and the delay between multicolor fields, 
IR+APT in particular, are conceptually equivalent, be- 
ing the result of interference between different multipho- 
ton pathways between the initial and final states. Ex- 
periments exercising this control are thus quite properly 
thought of as coherent control. 

This work is supported by the Chemical Sciences, Geo- 
science, and Biosciences Division, Office for Basic Energy 
Sciences, Office of Science, U.S. Department of Energy. 
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